Creation of multiple nanodots by single ions 
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In the challenging search for tools that are able to modify surfaces on the nanometer scale, 
heavy ions with energies of several 10 MeV are becoming more and more attractive. In contrast 
to slow ions where nuclear stopping is important and the energy is dissipated into a large volume 
in the crystal, in the high energy regime the stopping is due to electronic excitations only HJEHU. 
Because of the extremely local (< 1 nm) energy deposition with densities of up to 10 19 W/cm 2 , 
nanoscaled hillocks can be created under normal incidence 0J HJ |6). Usually, each nanodot is 
due to the impact of a single ion and the dots are randomly distributed. We demonstrate that 
multiple periodically spaced dots separated by a few 10 nanometers can be created by a single 
ion if the sample is irradiated under grazing angles of incidence. By varying this angle the 
number of dots can be controlled. 

As samples we used strontium titanate (SrTiOs) single crystals, a model substance for the class 
of perovskites (ABO3). The ferroelectric oxide SrTi03 is of use in a wide variety of applications: 
SrTi03 is used in catalysis and is discussed as a possible material for nuclear waste management 
ll8l . SrTi03 is the insulator of choice for the fabrication of future microelectronic devices. Because 
of its high dielectric constant thin films of SrTi03 have only a very low tunneling current compared 
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to conventionally used Si02 (UEI. In addition, it is widely used as substrate material for GaAs 
solar cells and for the deposition of high-temperature superconductors. For the latter, patterning of 
the substrate is known to influence the critical current densities IfTTTl . Therefore, the creation of regular 
nanometer sized structures in SrTiOs might as well proof useful for technological applications. 

After irradiation under grazing angles of a few degrees with respect to the surface plane (see fig. 0]) 
the samples were imaged by atomic force microscopy (AFM) in situ under ultra high vacuum (UHV) 
conditions. In fig. [T]we show an AFM image of a sample that was irradiated under two different 
incidence angles, resulting in two different types of defects. If irradiated at = 6°, slightly elongated 
hillocks occur. If the angle is changed to = 3°, round almost evenly spaced hillocks appear which 
are still overlapping. At angles smaller than = 3° the dots are clearly separated and the tracks look 
like pearls on a chain as can be seen in fig. [2] (this sample was irradiated at = l°and at = 2°). 
The number of chains corresponds to the number of ions. This means each chain containing dozens 
of nanodots is created by a single ion traveling through the solid. From analyzing more than 600 
individual chains irradiated under different angles we find that the length I of the chains can be 
controlled by varying the angle of incidence as can be seen in fig. [3] At = l°most of the chains 
are already about half a micron in length. The average height value we find for all chains stemming 
from irradiations at incidence angles of less than 3°is (2.5 ± 1.0) nm. They are of circular shape 
with a diameter (FWHM) of (25 ± 5) nm. The distance between individual hillocks depends on the 
angle of incidence and is e.g. (35 ± 3) nm for the marked chain in fig. [2] Similar measurements (not 
shown here) on irradiated Ti02 and AI2O3 show the same chains of dots as on SrTiOs. We take this 
as an indication that this phenomenon may be typical not only for perovskites but also for other oxidic 
materials. 

To date, the patterning of ferroelectric oxides has been achieved by self-organizing ('bottom-up') 
approaches lfT2l H3l as well as by 'top-down' methods such as lithography using low energy ions 
lfl4ll . evaporation through a shadow mask lfl5l or a combination of both lfl6l . With these techniques 
the created structures are either randomly distributed or shaped, or the periodicity is limited to 100 
nm at best. To our knowledge, no other technique that is able to create regular shaped nanodots on a 
perovskite surface with a spacing of 35 nm has so far been established. 

In order to describe the ion-solid interaction the so-called thermal spike (TS)-model has been 
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developed [171 . Toulemonde et al. applied it to the special case of swift heavy ions where electronic 
stopping dominates over nuclear stopping by orders of magnitude lfl"8l . It describes the formation 
of cylindrical tracks of amorphous or recrystallized material inside the solid. It does neither explain 
the formation mechanism of hillocks on the surface nor the formation of hillock chains. The starting 
point to gain an understanding of the formation of nanodot chains on SrTi03 are early experiments 
with van der Waals materials such as M0S2 or \VSe2. In these materials intermittent fission tracks 
were observed by transmission electron microscopy [19]. These discontinuous tracks were explained 
by the strong spatial anisotropy of a van der Waals crystal lattice. The n- electrons are fully contained 
within the crystal planes and there is only a weak hydrogen bond between the planes. Each time the 
projectile hits a lattice plane with a high density of electrons, enhanced electronic stopping occurs. 
This significant anisotropy leads to periodic spikes of radiation damage in this type of materials as 
observed by Morgan and Chadderton |fl9l . Such a Uiggered loss process could also explain the distinct 
features we observe. 

If the ion travels through SrTi03 it will continuously interact with the electronic system of the 
material but energy loss is more probable where the electron density is high. There, the sudden 
interaction of the projectile with the electronic system gives rise to a sharply peaked energy distribu- 
tion. How this electronic excitation finally leads to the modification of the material (craters, hillocks, 
tracks, etc.) is still not exactly known. The thermal spike model describes the energy transport out 
of the electronically heated region, whereas the Coulomb explosion model |[20l couples the electronic 
excitation to atomic motion caused by the repulsive forces acting in the transiently ionized region. 
Both models cover important aspects of the creation processes for material modifications Ell . In our 
experiment the strong ionization rate of the swift heavy ion leads to a significant charge imbalance. 
From Shimas empirical formula fl22l and extrapolating experimental data [23], respectively, we find 
that Xe ions moving through SrTi03 with 92 MeV have an effective equilibrium charge close to the 
original charge state of g g=23. However, our data does not allow us to make a clear correlation with 
neither of the two scenarios. In any case, for both models to be predictive the initial conditions need 
to be known. A strong anisotropy of the electronic structure as e.g. in a van der Waals material will 
influence the loss process significantly, be it connected either to a Coulomb explosion or to a thermal 
spike. 
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Note however, that SrTiOs is not a van der Waals material and the observed features cannot be 
simply attributed to the ion crossing crystal layers of a homogeneous electron density. If the anisotropy 
was exclusively parallel to the surface (as is the case with M0S2), we would observe that the number 
of nanodots was constant and the distance between dots would vary as a function of the angle of in- 
cidence. If the anisotropy was exclusively normal to the surface, we would observe a varying number 
of dots, but the distance between them would be constant. In our case the anisotropy is given parallel 
as well as normal to the surface. Therefore, the observed distance between dots and the number of 
the dots as a result of a triggered loss process needs to be discussed taking the full three dimensional 
electronic structure into account. To this end, we performed ab-initio density functional theory (DFT) 
calculations to determine the electronic density of SrTi03 and projected the density onto the plane of 
the traveling ion (see fig. [5]>. 

Most of the electrons are located around the oxygen atoms and the density is higher in the Ti02 
planes than in the SrO planes (see fig. @]). Simple geometrical considerations taking the electron 
density into account yield a very good qualitative agreement with the periodic defects occuring in 
SrTi03 as can be seen from fig. [5] In this approach we have used the DFT data to calculate the 
electron density encountered by an ion following a typical trajectory (azimuthal angle tp = 0°, see 
fig. |U) through a perfect SrTiOs crystal under a grazing angle of = 0.5°. On its way through a 
single lattice plane the ion interacts with the electrons of several dozens of oxygen atoms. Each peak 
in the fine structure (see upper panel in fig. [5]) is due to the interaction of the ion with the electrons 
surrounding one oxygen atom. The sum of these interactions within one crystal plane (large peaks in 
the lower panel of fig. [5]> creates a strong local excitation leading finally to a nanodot on the surface. 
The ion then typically travels several hundred Angstrom through the crystal without energy loss before 
it again encounters an area with a high enough electron density. 

Even though the exact mechanisms of the observed dot formation (fig. Q] and fig. 0) are not 
known, we may assume the space- and time dependent generation of excitation energy E s (f, t) along 
the ttajectory of the projectile r p {t) as well as the transport of excitation energy to play a key role 
in that process. In the energy regime considered here, the generation of excitation energy can be ap- 
proximately treated within the frame of the Lindhard model [24 ] of electronic stopping. Employing 
that framework, the excitation energy dE s {r p {t)) that is transferred from the kinetic energy of the 
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projectile into the electronic system within the time interval dt is proportional to the momentary ki- 
netic energy of the projectile and - even more important in our context - to the local electron density 
n e i(r p (t)) provided by the ab-initio DFT calculations as explained above. We assume the time evo- 
lution of the four-dimensional profile E(r, t) of excitation energy within the solid to be described by 
the diffusion equation 

dE dE 

— (f,t) - DV 2 E(f,t) = -^-(f p (t),E kin (f p (t)),n el (f p (t))) (1) 

with D denoting the diffusion coefficient. The projectile trajectory r p (t) entering Eq. (Q]) is obtained 
by numerically integrating the Newtonian equations of motion including an effective friction term 
originating from the Lindhard treatment. In view of the fact that on the one hand the lateral dimensions 
of the observed nanodots are of the order of several tens of nanometers, but, on the other hand, 
the electron density within one unit cell varies on the sub-Angstrom length-scale, a straight-forward 
numerical treatment of Eq. (Q]) is hampered by the complexity of the problem. 

It should be emphasized here that in our case (i) the consideration of the local electron density 
on a sub-Angstrom length scale and (ii) the necessity for the correct incorporation of the surface 
plane into the excitation energy transport process principally disallows the reduction of the problem 
to radial symmetry with respect to the swift heavy ion track as usually done in (TS)-model calculations 
ll25l . In contrast to the (TS)-model calculations, we propose to perform a nested two-step calculation 
approach as follows: In the first step, we take advantage of the straightness of the ion track and - 
for given impact parameters (6,0) - geometrically determine the set of intersection points through 
the unit cells along the trajectory of the projectile. This procedure is performed for a total trajectory 
length L which corresponds to a laterally projected range larger than the observed period length of the 
nanodots. Starting from the unit cell at the impact point, each traversed cell is discretized into more 
than 10 5 voxels to match the sub-Angstrom resolution of n e i(f,t). A numerical integration of the 
equations of motion of the projectile employing the highly-resolved electron density yields the total 
generated excitation energy dE s (r p (t)) as well as the traveling time dt for that cell. Naturally, the 
kinetic energy of the heavy ion is reduced after traversing one cell due to electronic stopping. 

This reduced energy as well as the exit point of one cell are taken as the new initial parameters for 
the analogous calculation carried out for the subsequent cell. Iteration of this algorithm yields a set 
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of (dE/dt)i which in the second step of our model, are treated as distinct point sources of excitation 
energy in a discrete representation of Eq. £T|) with a grid-spacing in units of one elementary cell. Using 
that coarse discretization the numerical solution of Eq. dTJ can be obtained via a Green's functions 
approach for the half space diffusion problem. 

In the present study we focus our interest on the striking periodicity of the nanodots which should 
already be present in the space dependence of the electronic stopping (ES) dE s /dx along the trajectory 
derived from ^§f-{r p {t)) of Eq. ([I]). Fig. [6] shows the result of an exemplary calculation of this ES 
for a 0.711 MeV/u Xe projectile hitting a SrTiOs crystal under = 0.5°and <fi = 10°. Two kinds 
of peaks can be seen, one originating from the Ti02 planes, the second one from the SrO planes. 
The average ES in the Ti02 planes is higher by a factor of 2 (integrating over the peak area). Thus, 
for this particular choice of impact parameters, the periodicity of the dot formation may be governed 
by the contribution of the Ti02 planes. The decrease of the maximum peak heights with increasing 
track length L due to electronic friction will be one factor limiting the total chain length /. A more 
quantitative discussion of features such as e.g. chain length, track radius and temporal dynamics goes 
beyond the scope of this paper because it would require a detailed analysis of the full four-dimensional 
excitation energy profile E(f, t) obtained from Eq. dD using the ES discussed above. 

Our model is further corroborated by plotting the measured length of the chains as a function of 
the angle of incidence as shown in fig. [3] The data can be fitted nicely by using /(©) = d/(tan@). 
Here, I is the length of the chain, and d is the maximum depth from where the excitation starts (see 
inset in fig. [3]). In addition, we frequently observe, that at the very end of a chain the height of the 
dots decreases monotonia These findings could be explained if we assume that the ion is traveling 
already too deep below the surface and the material modifications do not reach the surface any longer. 
In this way, we can determine the radius of the modified volume from our data. Our value of ca. 10 
nm is larger by a factor of 2 than the effective latent track radii derived from irradiation experiments 
on different insulators under 90° 12611271 . Repeating our experiment with an amorphous layer of Si02, 
we find periodic dots as well 1281 . It is not clear yet, what the origin of this periodicity is, but it could 
be the electron density, maybe even of the underlying crystalline Si substrate. 

In conclusion, we have demonstrated how to produce periodic nanodots on oxidic surfaces by a 
single ion hit. We propose that the anisotropic electron density of the material gives rise to a triggered 
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energy loss process. The resulting nanodots imaged by AFM thus represent a direct view of the pro- 
jection of the three-dimensional electronic density onto the surface. 
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Methods 



Experimental setup, sample preparation and image processing 

The experiments were performed at the beamline IRRSUD of the Grand Accelerateur National dTons Lourds 
(GANIL) in Caen, France. An UHV-AFM/STM (Omicron Nanotechnology, Taunusstein) was mounted directly 
to the chamber where the irradiation with swift heavy ions took place. The base pressure was 2 x 10~ 8 mbar 
in the differentially pumped irradiation chamber and 1 xlO~ 10 mbar in the AFM. As samples we used com- 
mercially available SrTiO3(100), TiO2(100) and Al2C»3(1102) crystals (CrysTec, Berlin). The samples were 
cleaned with acetone. We did not use any etching to get rid of excess SrO terminations since this would not 
play a role in the hillock formation. AFM images taken before irradiation to check the cleanliness of the samples 
frequently showed atomically fiat terraces separated by atomic steps. 

The samples were irradiated with a beam of Xe 23+ ions at 0.71 1 MeV/u scanned over the whole surface area 
of the target. Typically, a fiuence of 1 x 10 9 ions/cm 2 was chosen resulting in 10 impacts per /.im 2 on average. 
At this fiuence enough events are produced to obtain a good statistic but the probability that two ions hit the 
same spot resulting in non-linear phenomena is still sufficiently low. For SrTiOs the chosen projectile/energy 
combination results in an energy loss of 19.1 keV/nm as calculated with SRIM |29|. The angle of incidence 
with respect to the surface was varied between 0=l°and 0=6°. Due to the experimental setup the uncertainty of 
the absolute angle is not better than ±1°. Due to technical limitations the azimuthal angle <p was not controlled 
in our experiment. Immediately after irradiation all samples were measured with the AFM using the contact 
(F = 0.2 nN) as well as the dynamic (non-contact) mode (df = -10 Hz) in situ. Additional data from ex-situ 
measurements was used for fig. [3] No chemical etching or post-irradiation treatment was applied. This ensures 
that all observed topographical features can be attributed unambiguously to the irradiation. 

All AFM images were recorded with the Omicron SCALA software version 4. 1 and processed with the 
Nanotec Electronica SL WSxM software, version 4.0 Develop 8.3. From the raw data (400 x 400 data points) 
only a plane was subtracted. The colour code was changed using the palette flow.lut. No change of contrast (1) 
or brightness (0) was used. 

DFT calculations 

The density functional theory (DFT) calculations were performed using the ABINIT package 1 30 1 together with 
pseudopotentials generated by the fhi98pp code PD in order to determine the distribution of the electrons in 
SrTiC>3 . For the exchange-correlation energy the Perdew-Burke-Ernzerhof generalized-gradient approximation 
functional ||32l was used. In a first step, the equilibrium lattice constant of the investigated system was obtained. 
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The result is 7.53 ae for SrTi03. Subsequently the electron density at the equilibrium lattice constant was de- 
rived. A common kinetic-energy cutoff-energy of 96 Hartree for the expansion of the wave functions and a 
8x8x8 k-point mesh were used in all calculations. 
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Figure 1: Typical nanodots on SrTiOs created by irradiation with 92 MeV Xe ions. Frame size 500 x 
500 nm 2 , colour scale from to 4.4 nm. The sample was irradiated two times. The angle of incidence 
was = 3°and = 6°, respectively. Each dot or chain, respectively, corresponds to a single ion hit. 
To enhance the contrast false colouring was used. The arrow indicates the direction of the incoming 
ions. 
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Figure 2: Typical chains of nanodots on SrTiOs created by irradiation with 92 MeV Xe ions. Frame 
size 1300 x 1300 nm 2 , colour scale from to 3.7 nm. The sample was irradiated two times. The angle 
of incidence was = l°and 6 = 2°, respectively. The inset shows size and shape of the individual 
dots within a chain. 
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Figure 3: Measured length I of chains as a function of angle of incidence. The dash-dotted line is a fit 
according to our model with d = 9.8 nm (see text). 
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Figure 4: DFT calculation of the electron density (gray shading) of SrTiOs. Atom positions (besides 
the central titanium atom) are also visualized. The arrow indicates a possible projectile trajectory. 
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Figure 5: Simulation of the electron density which the ion encounters when traveling through the 
crystal lattice at a grazing angle of = 0.5° with respect to the (lOO)-surface plane of a SrTiOs single 
crystal (see fig. @]). 
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Figure 6: Electronic stopping dE s /dx along the track of a Xe projectile hitting a SrTi03 single crystal 
with 9 = 0.5°and = 10°. 
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